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ABSTRACT 

We have measured a velocity dispersion for the foreground galaxy in this 
gravitationally lensed system. Our dispersion confirms the prediction from 
lens models, provided that the source is distant enough {z^ > 0.2). Current 
interpretations of lensing statistics depend sensitively on how the optical and 
mass dispersions are related. For > 0.5, our observations favor (Jmass/cTopt — 1, 
but our uncertainties prevent us from clearly distinguishing between dynamical 
models. We could not obtain a unique mass profile from lensing constraints, but 
models with a constant mass-to-light ratio are possible. We also found unusual 
rotation in the lensing galaxy, and confirmed its redshift. 

Subject headings: Gravitational Lenses — Dark Matter — 
Radio Sources: individual (MG1549+3047) 
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1. Introduction 



Early examples of gravitational lensing (e.g., [Walsh etal. 1979|) consisted of several 



compact images of quasars. When the background source is extended (e.g., [Hewitt etal.\ 



1988[) , the images can join in a ring of emission around the lensing galaxy, and multiply- 
imaged source features can provide many constraints on the lensing mass distribution. For 
MG1549+3047, the background source is an extended radio lobe ( [Lehar etal. 1993]) , whose 



unusual structure is easily explained by a simple lensing model. Here we will be mostly 
concerned with the lensing galaxy, Gl (see Fig. 1 in [Lehar etal. 1993[ ), which lies directly 
in front of the Northwest radio lobe. 

Assuming an isothermal mass profile for Gl, Lehar etal. (1993) predicted its central 
velocity dispersion to be amass ~ 230 km s^^. A measured optical velocity dispersion aopt 
should be comparable to this value. Of course, neither the mass nor the stellar profiles 
need be isothermal. For a Hubble profile stellar population in an isothermal gravitational 



potential, Cmass/'^opt = \/3/2 ( [Gott 19T7[ ). This factor strongly affects the interpretation of 
lensing statistics (e.g., turner etal. 198^ ). However, considering a wider range of dynamical 



models, and taking observational limitations into account, Kochanek (1993) expects that 
o'mass/o'opt — 1 in most cases. Systems like MG1549-I-3047 provide a rare opportunity to 
test this prediction, since gravitational lensing can probe the total mass distribution. 



2. Observations 



We obtained long slit spectra of Gl on the night of 1992.04.28 at the 4.2 m William 
Herschel Telescope, using the red arm of the ISIS spectrograph. An EEV CCD detector 
was illuminated via the R600R grating, centered on 5750 A, with a 5300 A dichroic. This 
arrangement gave 0.73 A per pixel, or 38kms ^, and a spectral coverage of 5320 - 6180 A. 
The slit width was 1" on the sky (~ 3 pixels or 114kms~^), and was aligned with the 
major axis of Gl (PA = —151°). The field of view along the slit was ~ V. Six 30-minute 
exposures were taken. Using the same instrument configuration on 1992.05.01, we observed 
one radial velocity standard star of spectral type K2 (K2.1=HD155642), and two of type 
K5 (K5.1=HD155581, K5.2=HD156649). The spectra were reduced using standard IRAF 
procedures, and re-sampled into 44.6 km s~^ bins. Figure 1 shows the sum of all six Gl 
spectra (Gl.sum), along with a standard star spectrum. 

We measured the velocity dispersion in Gl using a standard cross-correlation method 
( [Ibnry fc Davis 1979| ; IRAF task "fxcor"). The three radial velocity standard star spectra 
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were used as templates. We removed 5th order spline fit baselines from the spectra and 
apodized by a 10% cosine bell, before Fourier transforming and applying a high-pass filter 
of ~ 10 in wavenumber. Finally, we fitted a Gaussian to the maximum correlation peak to 
determine its width and offset. We calibrated the correlation peak width using synthetic 
galaxy spectra (the template star spectra convolved with known Gaussian dispersions, 
see Figure 2). The results are summarized in Table 1. The K2.1 template gave higher 
cross-correlations than the K5 stars, so it is probably a better model of the Gl absorption 
features. We estimated the "random" contribution to the measurement uncertainty by 
performing separate cross-correlations for each of the six Gl exposures, and taking the 
rms of the dispersions. The "systematic" uncertainty, due to template mismatch, was 
estimated from the range of dispersions obtained when different stellar templates were used. 
Combining these estimates in quadrature, we found aopt = 227 ± 18kms~^. 

Our cross-correlation analysis also yields a radial velocity of cz = 33534 ± 103kms~^ 
for Gl. This result confirms a previously reported redshift of 0.111 ( Paturel etal. 19"89|) . 

Since Gl is extended, we extracted spectra for seven I'.'T apertures parallel to the 
galaxy trace in the co-added data. The cross-correlation results are given in Table 1, and 
Figure 3 shows the velocity profiles. The radial velocity shows a strong linear trend, which 
indicates that Gl is rotating. 



3. Discussion 



Our velocity dispersion, dopt = 227 ± 18kms~^, is very close to the (Tmass derived for an 
isothermal lens model ( [Lehar etal. 1993]) . The two mass estimates remain comparable for 
any source redshift > 0.2, confirming the basic lensing prediction. The radio properties 
and optical color of the source suggest that z^ > 0.3 (|Lehar etal. 19931) . Furthermore, 
MG1549-t-3047 is a steep spectrum source drawn from a flux-limited radio source list ( [Lehar 
1991| ). Assuming a spectral index of —1, the selection cutoff corresponds to the flux limit 
of the Parkes 2.7 GHz radio survey ( Dunlop fc Peacock 1990 ), for which {z) ~ 0.7 has been 
measured. Thus Zg > 0.5 is most likely. 

The systematic shift in radial velocity implies that Gl is rotating. The velocity is 
still increasing in the outermost apertures, so the maximum Vmax must exceed 200kms~^. 
Gl can be compared to an oblate isotropic rotator using the {vmax/o')* factor ([Kormendy 



|1982|) , which exceeds 2 . From its color and radial profile, Gl appears to be an early-type 
galaxy ( [Lehar etal. 19931) . Most ellipticals rotate more slowly, and this level of rotation is 
more characteristic of SO galaxies or the bulges in disk galaxies ( pavies etal. 198^ ). 
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The optical velocity dispersion can be compared with the total mass dispersion from 
lens models. Assuming a source redshift of = 1, the lens model gives an isothermal 
velocity dispersion amass = 230 ± llkms"^ ( [Lehar etal. 1995| ), which is the same as aopt- 



If the stellar tracers are assumed to have a Hubble profile, the expected optical dispersion 
will be 188 it 9kms~^ ( pott 19Y7|) ; and the strong rotation in Gl does not affect this. Our 
aopt differs from this prediction by almost two standard deviations, and agrees with the 
expectation of Kochanek (1993). However, (J mass incrCBiSGS by 25% if the source redshift 
Zg = 0.3 . So although our observation favors amass/ o^opt — 1, we cannot clearly distinguish 
between the two dynamical models. 

We can also probe the dark matter distribution by comparing the lens models with the 
optical surface brightness profile. Using the IRAF task "ellipse" , we extracted profiles from 
each of the CCD images reported in Lehar etal. (1993). Figure 4 shows the average profile 
for the / filter exposures. For sky surface brightnesses S(0), we used a model of the form: 

m = ^. [^ + {ifV , 

where P = 1 corresponds to an isothermal distribution, and P = 2 yields a Hubble profile. 
We found best fit models for each of the observed profiles, varying S^, 6c, and P. The best 
fit shape parameters, 6c = 0.24 ± 0.05 and P = 1.94 ± 0.05, did not vary with color. Here, 
the uncertainties are the rms of the parameter values, obtained from separate fits to each 
profile. For the mass distribution, we used a surface density model of the same form, with a 
quadrupole perturbation to account for the ellipticity e = (1 — 6/a) ( [Miralda-Escude 1993|) . 
The same image constraints and search method were applied as in Lehar etal. (1993). 
There is a strong correlation between 6c, e and P (see Figure 5). The values are very 
small, probably because the position uncertainties in Lehar etal. (1993) are overestimated. 
To estimate the parameter uncertainties, we performed Monte Carlo simulations, where the 
image positions were perturbed within their assumed O'.'l uncertainties. When P is close 
to its optical value, the best fit 6c is also consistent with the optical profile, provided that 
e ~ 0.3 . So although we cannot uniquely determine P, the lens constraints are consistent 
with the mass and light distributions being the same. 

We are now very close to a fully self-consistent lensing model for this system, and 
MG1549+3047 is promising for mass distribution studies. The background radio lobe has 
many structures, and it straddles the multiple-imaging boundary (see [Lehar etal. 1993| ). 



These properties are best exploited by the "Lensclean" method (^ocEaneF|^^ar^an 



19931) , which uses all of the source structure to constrain a lens model. The distance to 
the source is still unknown, however, and this deficiency affects many model predictions. 
Obtaining the source redshift remains an important objective. 
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Fig. 1. — Optical spectra for Gl (Gl.sum) and a template star HD 155642 (K2.1). The Gl 
spectrum is the sum from all six exposures. Wavelengths are given in Angstroms, and two 
major absorption features are labelled. 
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Fig. 2. — Velocity dispersion calibration curve. This shows the velocity dispersion 
(inkms~^) measured from model spectra with known "true" dispersions. The solid curve 
shows the result when K2.1 was used for both the model and template spectrum, and the 
dotted curves show other combinations of model and template. The heavy hne shows how 
the observed Gl correlation width converts to aopt- 
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Fig. 3. — Velocity profiles for Gl, as a function of aperture position (in arcsec NW of Gl). 
The lower panel shows the radial velocity (open circles) and the velocity dispersion (solid 
dots), both in kms^^. Error bars are the formal uncertainties from "fxcor". The upper 
panel shows the total counts in each aperture, and the extent of each aperture is shown as 
a horizontal bar. The dotted fines in both panels show the original 4'.'3 aperture. 
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Fig. 4. — Average surface brightness profiles for Gl, as a function of semimajor axis a (in 
arcsec), from the / filter CCD observations of Lehar etal. (1993). The lower panel shows 
the surface brightness E, in mag/arcsec~^. The middle and upper panels show the ellipticity 
e — (1 — b/a), and the major axis orientation (in degrees CCW from North). The error 
bars indicate the rms between profiles from each of the separate / filter exposures. Note that 
both e and 4> vary with radius, and that some of this may be due to atmospheric seeing. 
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Fig. 5. — Lens model parameters for various values of the mass index P. The lower panel 
shows the core radius 9c (in arcsec), with the best fit optical values indicated. The middle 
panel shows the ellipticity e — {1 — b/a), and the upper panel shows the P^r degree 
of freedom. The error bars show the 68% confidence intervals in 9c, from Monte-Carlo 
simulations (100 for each P). The dotted curves show the best fit parameter values if e < 0.2 
is applied as a constraint. 
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TABLE 1 

Spectral Cross-Correlation Results 



Spectrum Template Corr. r'^ cz^ <Jopt^ 



Summed Spectrum 



Gl.sum 


K2.1 


0.36 


17.1 


33534 


227 


Gl.sum 


K5.1 


0.34 


14.3 


33441 


215 


Gl.sum 


K5.2 


0.32 


13.9 


33435 


228 


Separate Spectra 










Gl.l 


K2.1 


0.25 


8.8 


33552 


200 


G1.2 


K2.1 


0.25 


9.4 


33527 


222 


G1.3 


K2.1 


0.22 


6.9 


33573 


213 


G1.4 


K2.1 


0.27 


9.0 


33537 


239 


G1.5 


K2.1 


0.25 


8.7 


33486 


221 


G1.6 


K2.1 


0.25 


8.9 


33533 


219 


Subapertures from 


Summed 


CCD 


Image 




+5'.'0NW 


K2.1 


0.08 


2.5 


(33806) 




+3'.'3 NW 


K2.1 


0.14 


4.5 


33831 


(109) 


+1'.'7NW 


K2.1 


0.22 


7.9 


33668 


207 


O'.'O NW 


K2.1 


0.35 


17.5 


33554 


226 


-I'.'TNW 


K2.1 


0.29 


9.9 


33450 


206 


-S'.'SNW 


K2.1 


0.15 


4.5 


33350 


(100) 


-S'.'ONW 


K2.1 


0.06 


1.8 


(33432) 



^Tonry &l Davis (1979) signal/noise estimator. 
^Velocity in (kms~^) 



